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Abstract. Brightwell, van den Heuvel and Stougie proved that the diameter 
of an m X n transportation polytope is at most 8(m + n — 2), a factor of eight 
away from the Hirsch Conjecture. This bound was improved to 3(m + n — 1) 
by Hurkens. We investigate diameters for certain classes of transportation 
polytopes. 



1. Introduction 

In |SanlO| . Santos provided the first counterexample to the famous Hirsch Con- 
jecture, which asserted that the diameter of the 1-skeleton of a <5-dimensional poly- 
tope with facets is at most <f> — 5. This bound, which we call the Hirsch bound, 
was however natural and plausible. Although there is little hope to characterize 
the polytopes which satisfy this bound, it is interesting to investigate families of 
polytopes for which the bound may hold. The transportation polytopes constitute 
one candidate for such a family. 

Consider a supply function m : M —> M>o on a set Mofra sources and a demand 
function n : N — >• M>o on a set N of n sinks. The m x n transportation polytope 
-P m> n is a convex polytope obtained by intersecting the positive orthant M^f xAr with 
the following m + n afhne hyperplanes: 

We always assume that Xl At m (/ i ) = Si/ n ( I/ )j such that P m ,n is non-empty. Intu- 
itively, a point of the transportation polytope P m . n is an assignment of quantities 
to be transported between the sources of M and the sinks of N on each edge of the 
complete bipartite graph Km,n such that the total quantity that a source \i pro- 
vides corresponds to its supply m(fi) while the total quantity that a sink v receives 
corresponds to its demand n(^). Optimizing transportation costs naturally gives 
rise to linear optimization problems on transportation polytopes, and thus leads to 
the question to evaluate the diameter of transportation polytopes. 

Any non-empty m x n transportation polytope has dimension (m — l)(n — 1) 
and at most mn many facets. Consequently, if the Hirsch Conjecture is true for 
transportation polytopes, the diameter of any m x n transportation polytope will 
be at most (j>— 6<m + n— 1. In BvdHS06 Brightwell, van den Heuvel, and 
Stougie gave a bound of 8(m + n — 2), which was then improved to 3(m + n — 1) by 
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Hurkens [Hur . However, it is not clear so far whether the transportation polytopes 
satisfy the Hirsch bound or not. 

In this note, we study the diameter of a specific subfamily of transportation 
polytopes. Given a degree function d : M N>o, we consider the m x n trans- 
portation polytope Pd : = -P m ,n where the supply and demand functions are defined 
by m(/i) := l + md(/x) and n(z/) := m. (It forces n = 1+X)u to ensure that to- 
tal supply equals total demand.) We show that the diameter of Pd does not exceed 
twice the Hirsch bound. Furthermore, when d = 1, the polytope Pd is an mx (1+m) 
generalized Birkhoff polytope, and we prove that its diameter is precisely given by 
the Hirsch bound. Along with the diameters of Birkhoff polytopes [YKK84J and a 
recent extension to partition polytopes by Borgwardt Borlf, this provides one of 
the first subfamilies of transportation polytopes with this property. 

Our study of the polytope Pd essentially relies on its interesting combinatorial 
structure: the 1-skeleton of Pd is (isomorphic to) the pivoting graph on spanning 
trees of the complete bipartite graph Km.n where each source /i has prescribed 
degree 1 + d(/i). We first study independently this family of spanning trees, and 
relate it later on with the transportation polytope Pd. 

2. Bipartite spanning trees with specified source degrees 

Denote by Km,n the complete bipartite graph between a set M of m sources and 
a set N of n sinks. Fix a function d : M — > N>o and assume that n = 1 + ^2 d(/i). 
We consider the set Td of spanning trees of Km.n where every source [i has de- 
gree 1 + d(/i). Note that we prescribe the source degrees of these trees, but not 
their sink degrees. Any such tree has e := 1 + d(/i) edges. The value of n was 
chosen so that e = m + n — 1, and thus Td is nonempty. We first give a closed 
formula for the cardinality of Td- 

Proposition 2.1. The number of trees in Td is precisely 

(n-iy.n™- 1 
' dl n,d(/x)! ' 

Proof. We adapt Andre Joyal's proof of Cayley's formula. We fix a total order -< 
on M and a special sink v* E N. Let IZd be the set of trees of Td oriented towards 
sink and with a marked sink (which may or may not coincide with v„). The 
number of trees in IZd is n|7d|- Let Sd be the set of directed subgraphs of Km,n 
where every source u. has in-degree d(/x), and every source and sink has out-degree 1, 
except j/*, which has out-degree 0. Since X^d(^) = n — 1, the number of graphs 

in Sd is TpaTOl • ^ e now define inverse bijections $ : Sd — > T^d an d * : TZd $d- 
Consider a directed graph S in Sd- The connected component Si of S containing 
sink is a tree oriented towards . The other connected components S2 , • ■ ■ , S p 
of S are graphs of functions: each is formed by a directed (source-sink alternating) 
cycle together with some rooted trees which are oriented towards their roots and 
glued to the cycle by their roots. For k > 2, we denote by fik the source of 
the cycle of Sk which is maximal for and by the sink following fi^ in this 
cycle (fix also v\ = v*). By reordering, we can assume without loss of generality 
that fj,2 ~< ■■ ■ ~< (J-p- See Figure [TJ To obtain the oriented tree <&(S) from S, we 
open the cycle in each by erasing the arc between u. k and v>k, we concatenate 
the resulting trees by adding an arc between and and we mark sink v p . In 
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other words, 

:= (S \ {( Pk , u k )\ke{2,... ,p}}) U {(wt, | k G {2, . . . ,p}} 
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Si 

Figure 1. Schematic decomposition of S G Sd into Si,...,S p . 
The map $ removes the long blue arcs and inserts the dashed green 
arcs. 



Consider an oriented tree R in TZ A with marked sink v* (which may or may 
not be v^). Consider the path from to v* (note that this path goes against the 
orientation on R). Denote by the source succeeding v*, and inductively choose 
/is, . . . , fx p so that Hi is for all i the first source on the path that is greater than 
with respect to -<. For k < p — 1, let v k be the sink preceding /Zfc+i along this path, 
and define i/ p := z/\ See Figure [2j To obtain the directed graph ^(R) from R, 
we disconnect R by erasing the arc between \i k and v k _\ and we create cycles by 
adding an arc between fi k and v k . In other words, 

*(R) ■= (R\{(p ki u k -!) | k G {2, . . . ,p}}) U {(m*, i/*) | fe e {2, . . . . 
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Figure 2. Schematic depiction of i? G TvLd with marked sink i/*. 
The map $ removes the long blue arcs and inserts the dashed green 
arcs. 



It is clear that $ and \& are inverse bijections. We thus obtain that 

\K A \ \S A \ {n-\)\n m - x 



n u d(/i)! 



□ 



When d = d G N>o is constant on M, the number of trees in 7d is — T^ya — , 
which recovers a formula of Klee and Witzgall (see |KW68j or Theorem 2.4 in 
Chapter 6 of [YKK84 ). This counting can alternatively be derived from classi- 
cal generatingfunctionology. For that, consider the collection of all bipartite trees 
whose sources are unlabeled and have degree d, whose sinks are labeled, and which 
are rooted at an arbitrary sink. Denote by Od{n) the number of such trees with n 
sinks, and let Qd(z) ■— J2 n ®d{n)z n /n\ be the corresponding exponential generating 
function. According to the dictionary of generating functions [FS09] , this function 
satisfies Qd{z) = z(d(®d(z)), where Q(y) — exp(y d /d!). The Lagrange inversion 
formula ensures that 0d{n) = n\ [z n }<d d (z) = (n — 1)! [y n ~ 1 ]C,d{y) n , where [x p ]f(x) 
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denotes the coefficient of x p in the series expansion of the analytic function f(x). 
Consequently, the number of trees of 7d when d = d is given by 

mW d (l + md) = m\(md)\ md 1+md = (md)!(l + md)^ 1 = {n - l)\n m ^ 
l + md l + md W MW (dl) m (d\) m 

In fact, this method based on generating functions can also be extended to obtain 
the formula of Lemma |2.1| for general d, using the multivariate version of Lagrange 
inversion formula [BR98 . 

Let T be a tree in 7d and e be an edge of Km,n not in T . There is a unique cycle 
in TU{e}. Let e' be the other edge in this cycle adjacent to the source of e. Then T 
and T' ■= (T \ {e}) U {e'} are the only two trees in 7d contained in T U {e}. We say 
that T" is obtained from T by pivoting on e. Note that this operation removes the 
edge e' from T. 

We consider the pivoting graph whose vertices are the trees in 7d, and whose 
edges are pairs of trees which differ by pivoting on an edge. We study the diameter 

of Gd- 

Theorem 2.2. The diameter ofGd is at most 2J^ d(/z) = 2n — 2. 

Proof. Consider two trees S and T in 7d- We say that a source is consistent in (S, T) 
when all edges incident to it coincide in S and T . We will prove that we can make 
all sources consistent after at most 2 d(/i) pivots in S and T, which means that 
the distance between S and T in is at most 2 We proceed by induction 

on the number of non-consistent sources. 

Assume that S and T are distinct so that there exists at least one non-consistent 
source in (S, T) . We pick any source /i* and choose a non-consistent source \x with 
maximal distance to fj,* in S. Removing fi from S results in l+d(/i) components. We 
denote by Vq the vertex set of the component containing /i* and by v\, . . . , V^aO 
the vertex sets of the other components. We also denote by Vk the sink in Vk 
adjacent to /i in S. The choice of /i ensures that all sources among V\, ... , Var^) are 
consistent in (S, T) . In particular V\ , . . . , Va^) induce connected subgraphs of T. 

In at most 2d(/x) pivots, we will transform S into S and T into T such that \i 
becomes consistent in (S, T) while ensuring that sources that are consistent in (S, T) 
stay consistent in (S,T). 

We will perform our pivots in two phases. For the first phase, let us assume 
that there is an edge (fi,fk) of <S* that is not shared by T. The unique cycle in 
TU {(/U, fk)} contains a second edge (fi, v) incident with \i. We distinguish two 
cases. If v ^ {vi, . . . , Vd(fi)} then we perform a pivot on (p, Vk) in T, which removes 
(/i, v) from T. 

In the other case, when v — vi for some £ € {1, . . . , d(/i)}, we observe that 
the unique cycle in T U {(/z, Vk)} contains a non-consistent source a other than /x; 
otherwise the cycle would be contained in the tree S. Next, the source /i has a 
neighbour r ^ {vi, . . . , t'dO)} m T. See Figure |3| We now perform first a pivot on 
(a, t) in T. In the tree T" obtained from T by this pivot, we find that inserting the 
edge (/ijffc) results in a cycle that contains (h,t). Consequently, our next pivot, 
which we perform on (u.,Vk) in T', removes (/x,r). We continue with this process 
until all the edges (fi, Vk) lie in T. This concludes the first phase. 

After the first phase, either [i has become consistent or there is exactly one edge 
left in T that is not shared by S. Clearly, the source vertex of this edge is fi, while 
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Figure 3. Phase one pivots on T: the cycle in T U (/it, v^) is de- 
picted in red. Pivoting on (cr, r) in T removes the edge to the right 
of <7, producing T' . Then pivoting on (/x, Vf.) m T' removes (/x,t). 



its sink, which we denote by k, is contained in Vq. Then pivoting in S on (/x, k) 
makes \x consistent. 

The previous arguments show that the source /x can be made consistent in at 
most 2d(/x) + 1 pivots: we need at most two pivots to insert (/x,!^) in T for 
each fc > 1, and one final pivot to insert (/x, k) in S. To obtain the claimed re- 
sult, we observe that the current procedure requires 2d(/x) + 1 pivots only in the 
case where we start with no edge between /x and v\, . . . m T. But in this 

case, we can drop the last 1 since at the very beginning, we only need one pivot to 
insert (/x, V\) in S. □ 

3. Background on transportation polytopes 

Given a supply function m : M — > K>o on a set M of m sources and a demand 
function n : N — > M>o on a set N of n sinks, the transportation polytopc P m rl has 
the following inequality description: 

(1) V/x, Vi/, af^.i, > 0, ^2 v , x^y = m(/x), and ^ cc^^ = n(i/). 

Again, we assume that P m) n is non-empty, i.e. that m(/x) = X)i/ n (^)- Since 
P m n is defined by m + rx — 1 linearly independent equations in an mn-dimensional 
ambient space, it has dimension 5 ■— (m — l)(n — 1). 

From the inequality description ([!]), the mxn transportation polytope P mi n has 
at most "ran facets, each defined by an inequality of the form x^ u > 0. Reciprocally, 
the following lemma affirms that an inequality x„ v defines a facet as soon as it is 
tight. 

Lemma 3.1. Let mn > 4. The inequality x^ v > defines a facet o/P mi „ if and 
only if there is an x £ P m ,n such that x^. v = 0. 

This follows from e.g. Theorem 3.1 in Chapter 6 of |YKK84j . which states: 

Proposition 3.2 f |YKK84) . pg. 271). Let mn > 4. The inequality x^*^* > 
defines a facet of P m , n if and only if m(/x*) + n(^*) < m (A*) 
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Let x G P m ,n- The support of x is the subgraph supp(x) of Km,n consisting 
of the edges (/x, v) for which x^ M > 0. The point x lies in mn — | supp(ir)| many 
facets of P m ,n- In this article, we assume that P mn is simple, i.e. that every 
vertex has precisely m + n — l non-zero coordinates, or equivalently, that no proper 
subsets M' C M and N' C N satisfy EueM' m (/ i ) = TlveN' n ( v )- The following 
proposition summarizes the properties of the supports of the vertices of P mi n (for 
a proof, see e.g. |YKK84j ). 

Proposition 3.3. Let P mn be a simple transportation polytope. 
(i) A point x G P m .n is a vertex of the polytope P m ,n if and only i/supp(x) is a 

spanning tree of Km,n- 
(ii) A vertex x o/P mn is determined by its support supp(x). 
(in) The supports of two adjacent vertices of Pa differ in precisely two edges. 

We close this background section with a relevant family of transportation poly- 
topes. The Birkhoff polytope B m of size m is the transportation polytope whose 
supply and demand functions are both constant to m. Its vertices are precisely the 
permutation matrices. Note that the support of a permutation matrix is a perfect 
matching, so that the classical Birkhoff polytope is not simple. We also want to 
underline that the Birkhoff polytope B m is known to have diameter exactly 2 (see 
Theorem 1.7 of Chapter 5 in [YKK84 ). The m x n generalized Birkhoff poly- 
tope B miTl is the transportation polytope whose supply function is constant to n 
and whose demand function is constant to m. Note that it is a simple polytope if 
and only if m and n are relatively prime. As a consequence of the next section, we 
will obtain the diameter of generalized Birkhoff polytopes B m i +mc i for m,d> 1. 

4. A SPECIAL CLASS OF TRANSPORTATION POLYTOPES 

In this section, we study a specific family of transportation polytopes whose 
combinatorial properties are closely related to that of the family Td of trees defined 
in Section [21 We continue to use the notation of Section [21 

Let m : M — > M>o be the supply function defined by m(/i) = 1 + md(/i) for 
all /i <E M, and let n : N — > E>o be the demand function defined by n(^) = m 
for all v G N. Then the m x n transportation polytope Pd ■— P m , n is non-empty 
(since m (^) = n (^)) and simple (since E M eM' m (M) ¥= Er/eJV' n ( v ) for an y 
proper subsets M' C M and N' C N). 

Example 1. If d = d, then n = 1 + md and we have m = n and n = m. The 
polytope Pd is thus the generalized Birkhoff polytope B m n of size m x n. Note 
that m and n = 1 + md are relatively prime, so that Pd = P m ,n is simple. 

The following lemma relates the family Td with the transportation polytope Pj. 

Lemma 4.1. The support function x i— > supp(x) is a bijection from the vertex set 
of Pd to the set Td of trees, which sends the 1-skeleton of Pd to the graph Qd- 

Proof. Let T = supp(x) for a vertex x of Pd- The degree of a source /i in T is at 
least d(//) for any source /i. Since T has precisely m + n — l = E M edges, each 
source fi has degree exactly d(/i) in T, and thus T is in Td- 

We prove the reciprocal statement by induction on m + n. Let T € 7d- We need 
to exhibit a point x G P m ,n whose support is T. Note that x will automatically be 
a vertex of P m! „ since its support is a tree (see e.g. |YKK84| V 
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Assume first that a sink r is a leaf in T only adjacent to a source a. Let N' = N \ {t}, 
let T be the restriction of T to MUN', and define d' : M N> by d'(cr) = d(<r)-l 
and d'(/x) = d(/i) if /i 7^ a. Since T' G 7d', the induction hypothesis ensures 
the existence of a vertex x' G C M. M * N whose support is T' . We define a 
point x G $i MxN by extending x' with i ffjT = n(r) = to and x^ hT — for \l ^ a. 
Then x is a vertex of Pd with support T. 

Assume now that a source a is a leaf in T only adjacent to a sink r. Let M ' = M \ {a}, 
let T' be the restriction of T to M' U iV, and let d' be the restriction of d 
to M' . Since T' G 7d', the induction hypothesis ensures the existence of a ver- 
tex x' G -Pd' C M xjv whose support is T'. We orient the edges of T' towards r. 
We define a point a; G R MxN by: 

• = a;' v if there is an arc oriented from ji to v in T", 

• 2^,1/ = !/ + 1 if there is an arc oriented from v to /i in T", 

• = if (/i, ^) is not an edge of T. 

Then a; is a vertex of Pd with support T. 

We thus proved that the support function x i— > supp(x) is a bijection from the 
vertices of Pd to 7d- For the second part of the statement, observe that two vertices 
of Pd are adjacent if and only if their supports differ in precisely two edges, i.e. in 
a pivot. □ 

More generally, the support function provides an isomorphism from the face 
lattice of the polar polytope of Pd to the simplicial complex whose maximal faces 
are the complements of the trees in 7d- 

Lemma 4.2. // at least two sources have prescribed degree at least 2 (in other 
words, if |d _1 (N>i)| > 2), then the number 4> of facets of Pd is precisely mn. 

Proof. By Lemma [3. 1[ the inequality x^^ > defines a facet of Pd if and only if 
there exists a vertex x G Pd with x^^ — 0. Such a vertex exists because there is a 
tree T G 7d with v) £ T since 1 + d(p) < n. □ 

In the conditions of this lemma, the Hirsch bound for the diameter is thus 
H ■= cf) — 5 = to + n — 1. Applying the result of Theorem |2.2[ we obtain: 

Theorem 4.3. The diameter of Pd is at most 2n — 2 < 2H . In particular, if d = 1, 
then the diameter of Pd = P mj i+ m is at most the Hirsch bound H. 

Note that the generalized Birkhoff polytope P 2 ,3 is an hexagon. Consequently, 
its diameter is strictly less than the Hirsch bound H . This situation does not occur 
for larger to x (1 + to) generalized Birkhoff polytopes: 

Corollary 4.4. For any to > 3, the diameter of the tox (1+to) generalized Birkhoff 
polytope B m ^i +m coincides with the Hirsch bound 2m. 

Proof. The upper bound is given by the previous corollary. For the lower bound, 
it is easy to construct two trees of 7d with disjoint edge sets. □ 

Pak [PakQO] investigates the /-vector for to x n generalized Birkhoff polytopes 
when n = m + 1. Using Proposition |2.1| and the simplicity of Pd, we obtain the 
first two values of the /-vector for a slightly more general class of transportation 
polytopes: 
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Corollary 4.5. The number of vertices and edges of Pd are respectively 

/o(Pd) ~ n M d(M)! and h[Pd) - 2n^p ' 
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